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Abstract 

(N , 

. We show that the Painleve equations Pm — Pyi can be derived (in a unified way) from a pe- 

riodic sequence of Darboux transformations for a Schrodinger problem with quadratic eigenvalue 
\ dependency. The general problem naturally divides into three different branches, each described 
' by an infinite chain of equations. The Painleve equations are obtained by closing the chain peri- 
■ odically at the lowest nontrivial level(s). The chains provide "symmetric forms" for the Painleve 



equations, from which Hirota bilinear forms and Lax pairs are derived. In this paper (Part 1) we 
analyze in detail the cases Pm — Py, while Pyi will be studied in Part 2. 

3 ■ 

^ : 1 Introduction 



The aim of this paper is to provide a bridge between two extremely fruitful descriptions of 
the Painleve equations in the context of integrable systems: the (Darboux) dressing chain 
approach of the Russian school [H 12 El IH El E] and the tau-function based approach of the 
Japanese school [7lElEllIDllIIlliaiIllIlllI3lini- 

In the seminal paper j3] Adler proposed a complete description of the Painleve equa- 
tions Pjj-vi by means of so-called Darboux chains. These generalize the usual dressing 
chain |1| which, in its periodic case, was already known to provide a description of the 
Piv and Py equations Adler succeeded in describing the Schlesinger-Backlund trans- 
formations and notably the Weyl-group action of these transformations for the Painleve 
equations, merely by starting from the Darboux transformations that underlie their asso- 
ciated (periodic) chain equations. However, its relationship with Okamoto's theory of the 
Painleve equations (in which Weyl-groups play a central role) has still not been elaborated. 
In the Okamoto description, the Weyl-group action of the Backlund transformations for 
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the Painleve equations is closely related to the notion of a tau-function jH 13 UHl UH 1121 UH] 

(see for an overview), which in turn is connected to the existence of isomonodromy 
problems associated with the Painleve equations [H]. In recent years this (afiine) Weyl- 
group structure of certain Backlund transformations for the Painleve equations, together 
with the tau-function description of the Hamiltonians underlying these equations (see JHl 
for a survey of results), has given rise to a "symmetry-based" approach to the Painleve 
equations and related systems [HI Hn] . In [Hj Noumi and Yamada proposed a systematic 
description of a class of dynamical systems comprising the Pjv and Py equations, which 
all possess Backlund transformations that make up an affine Weyl-group and which allow 
for tau-function descriptions of their Hamiltonian structures. So called "symmetric forms" 
of these dynamical systems play a central role in this approach |16j . 

We will show, using the associated Hirota bilinear representations, that periodic closings 
of the standard dressing chain considered in [21 El E] ^ and from which Pjv-v are obtained 
at periods 3 and 4 respectively - correspond exactly to the dynamical systems introduced 
in Our approach is, however, not restricted to the Piv-v equations or similar cases. 
We shall also explicitly derive a symmetric form of the Pm equation and of a whole class of 
related equations (which first appeared in [1|), all of which can be cast into Hirota bilinear 
form, i.e., their solutions can always be expressed in terms of tau-functions. 

A major part of this paper is devoted to a systematization of the results obtained 
by Adler [3] in order to be able to derive all relevant Darboux (dressing) chains from a 
common starting point and to analyze their properties using the same techniques. This 
unified approach then not only pays off when bilinearizing the Darboux chains we obtain, 
but also and especially so when deriving Lax representations of those chains. Although a 
lot of attention has been devoted to the spectral properties of the linear equations whose 
Darboux transformations underlie the various (Darboux) dressing chains associated to the 
Painleve equations 2112011^1; no Lax representations of these chains have been offered in the 
literature pertaining to Darboux chains, although the Noumi- Yamada symmetry approach 
does yield such representations. In the following we shall, systematically, derive Lax pairs 
for the various Darboux chains presented throughout this paper. In the Piv-v case these 
Lax pairs will turn out to be identical to those obtained in [TBI 1^ , but we believe the 
result concerning Pjjj is new. 

The structure of the paper is as follows: We start by defining the Darboux trans- 
formation for a Schrodinger operator with an energy-dependent potential and derive its 
associated chain equation. The analysis naturally divides into three branches covering 
Pjv — Py, Piii, and Pyi, respectively. In Section 3 we study the Pjy — Py branch in 
general and then in Sec. 4 give the details for Pjv and in Sec. 5 for Py. Then we study the 
Pjii branch in general and in Sec. 7 give the specific details for Pm itself. As mentioned 
above, the results concerning the Pyi case will be discussed in Part 2. 

While finishing this paper we came across Ref. [3] by Shabat, which is also concerned 
with systematizing the Darboux-Painleve connection. It is therefore not surprising that 
our analysis runs largely parallel with that paper and that many formulae are identical or 
similar. However, as we wish to focus on the specific connections to the Painleve equations, 
continuing with their bilinearization and Lax pairs, our presentation is both more specific 
and more detailed. We believe that in doing so the method and its internal machinery gain 
further clarity. 
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2 The Darboux transformation 



Roughly speaking, a Darboux transformation is a transformation between two hnear (or- 
dinary or partial) differential equations, of the same overall appearance but for different 
values of their coefficients: the solutions of the first equation will be mapped to the solu- 
tions of the second, and the changes in the coefficients of the equation are also provided 
by the transformation (see, e.g., [22] for a classic introduction in the context of integrable 
systems). We are interested in chains of such transformations applied to the Schrodinger 
equation. 

2.1 The general form of the chain equations 

Consider the spectral problem associated with the Schrodinger operator ^ 

Lj{u, \)tpj{\, x) = 0, where Lj{u, \) := + Uj{\, x). (1) 

Here j indexes the eigenproblem, and A are the eigenvalues (assumed non-degenerate). 
Note that we make no assumptions about boundedness or normalizability of the eigenfunc- 
tions. We shall refer to the functions Uj{X,x) as "generalized potentials" as they include 
dependence on the eigenvalues A. Next define the operator 

G,(A,x) :=A,(A,x)(9,-F,(A,x)), (2) 

and then using this define new functions 

V'j+i(A,x) := Gj{X,x)ipj{X,x), (3) 

for each eigenfunction ipj{X, x) of the original problem at level j. The new functions ipj+i 
will then be eigenfunctions of a new operator Lj^i (with the same, generic eigenvalue A), 
provided that L^+i satisfies the operator identity 

Lj^iGj{\,x) = Gj{\,x)Lj{u, \), (4) 

for some Gj{X,x) := Aj{X,x){dx — Fj{\,x)). Assuming that Lj^i is a differential operator 
(i.e., polynomial in d^) one finds that 

F,(A, x) = F,(A, x) - 2(log A,(A, x))', (5) 

Lj+i = Lj{uj+i,X), (6) 

Uj+i{X, x) = Uj{X, x) + [2Fj{X, x)Aj{X, x) - Aj{X, a;)']'Mj(A, x), (7) 

Fj{X, x)' + Fj{X, xf + UjiX, x) = fij{X)Aj{X, x)-^, (8) 

where fj.j is an integration constant (the '-notation stands for ^). 

The transformation (jH)) from ipj to ipj+i will be called a Darboux transforma- 
tion iff at each step j the operator Gj is such that it annihilates some chosen 
eigenfunction (pj of (HJ having eigenvalue Uj. 
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This necessarily implies that Fjivj^x) = {logipj)x and therefore these Fj also satisfy the 
equation 

Fj{uj,xy + Fj{uj, x f + Uj{uj, x) = 0. (9) 

This should be compatible with (jS)) and therefore fijii^j) = 0. 

Now subtracting (jSI) for j and for j ' + 1 and using ((7j) to eliminate the potentials u we 
get a chain of equations 

F,+i(A,x)' + (F,(A,x) - {\ogA,{X,x))'y+ 

F^+i{X,xf - {F,{X,x) - {\ogA,{X,x)rf+ 

/i,(A)A,(A,x)-2-/i,.+i(A)A,+i(A,x)-2 = 0. (10) 

In the following we shall think of this chain equation as a generator for A-independent chain 
equations (the dressing chains). Thus the main equations underlying our analysis will be 
((Tj), which gives the change in u once F, A are given, and (fTU)). the chain equation. 

2.2 Linear problem for the chain equations 

The commutation relation (@|) which generates the chain equation (jlUj) can also be rewritten 
as 

- = or A'^L,+,G, = G.A^L,, (11) 

which allows one to interpret the system ()1|3|) as a kind of linear problem for the chain 
equation. Let us define 

Mj := dx - AJ^S - Fj, where S : tpjiX, x) ^-^ ^j+i{X, x), (12) 

then since 

Mjtljj{X,x) = (13) 

(i.e., equation (jS))) we can eliminate from the eigenvalue problem (PJ and obtain the 
following second order difference equation for the eigenf unctions iljj{X,x): 

L^^,{X,x)=0, (14) 

where 

:= AJ^S Aj'S + Aj'S F, + Fj Aj'S - (log Aj)' Aj^S + fi.Af (15) 

(see also [11123 1211 )• The compatibility condition of the system p3ll4p . which is of course 
a natural consequence of the compatibility of Eqs. ((T)) and ©, takes the form: 

A, {M,L^ - LjM,) + 2A^Lj = 0, (16) 

and is satisfied iff ()10|) holds. This linear problem will play an important role later on when 
we derive explicit forms for the Lax pairs for various reductions of the dressing chains. 
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2.3 The specialization that contains Pjjj — Pyi 

The above is as far as we will go with the general setting. Now we specialize to the 
A-dependence 0j 

Uj{X, x) = — A^ + Xvj{x) + Wj{x) (17) 

for the potential, and 

F^{\,x) = \hj{x) + f,{x) (18) 

for the Darboux transformation. We will now show that this simple case is still rich enough 
to allow for three different infinite hierarchies of equations, containing Painleve equations 
as their simplest nontrivial members. 

With the above assumptions we find from (0) and (jH)) that both (logy4j(A, x))' and 
Aj{X, x)~^ should be polynomials in A. This is only possible if the A-dependence in Aj{X, x) 
is multiplicative, in which case it can always be incorporated into the integration constants 
/Xj(A) in (jSJ. Furthermore, comparing (jHJ and ^ shows that should contain an 

overall factor of A — uj and thus we can take it to have the form 

/i(j,A) = (A-z.,)(Aa,+/5,). (19) 

When the above ansatze are substituted into equation (jSJ we find, at successive powers 
of A, the following equations: 

A^ : Ajixy^aj = hj{xf - 1, (20) 

A : A,{x)-\P, - a,u,) = h.ix)' + 2f,{x)h,{x) + v,{x), (21) 

1 : -Aj(x)-2/5jZ/j- = f){x)' + fj{xf + Wj{x). (22) 

The solutions of this set of equations naturally split into 3 distinct branches: 

1. If /i^ = 1 then we can take (without loss of generality) hj = 1, <yj=0 and /3j = 2, and 
solve Aj'^ = fj{x) + ^Vj{x) from (j2T}. The chain equations (|Tn|l obtained from this 
branch will contain Pjjj ClS db special case. 

2. If hj = we can take aj = —l,Aj{x) = 1 and equation (|2ip then becomes Vj{x) = 
Uj + Pj. For convenience we set Vj{x) = and Pj = —Uj, which implies that /i(j. A) = 

— A^ and Uj{X,x) = — A^ + Wj and we therefore obtain the usual dressing chain 
for the Schrodinger equation fSlE], after changing notation A^ — > A, z/| — > uj. This 
branch contains Pjv and Py- 

3. Finally, in the generic case we choose aj = —1 so as to get Aj{x)~^ = 1 — hj^x)"^. 
Equations ()21|) and (j^^ then yield ODEs for h and /. In Part 2, this branch will be 
shown to contain Pyi- 

Note also that the Pjv,Pv branch is obtained as a limit /i — ^ from the Pyi branch. The 
Pjjj branch can also be obtained as a limit from the Pyj branch, but the limit is singular: 
Let 

hj = 1 - 6{fj + ^Vj), aj = -2e, (23) 

and then as e — > the leading terms from equations ()2UI21|) yield Aj"^ = fj + ^vj and 
Pj = 2, respectively. 

In all these branches the Painleve equations are obtained when the derived chain equa- 
tions are closed (with "more or less" periodic boundary conditions as we shall see later 
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on) after 2, 3 or 4 steps. Each time there exists a "hierarchy" of equations, obtained by 
closing the chain after a higher number of steps, but we will not discuss these higher order 
equations in the present paper. 

3 The Pjv-V branch in general 
3.1 The generic chain equations 

The simplest realization of the above scheme is found in the case of the (ordinary) Schrodinger 
equation {dl+w{x) — \)ip{\,x) = 0, i.e., for a generalized potential ^^(A,^) in the operator 
© of the type 

Uj{X,x) = Wj{x) - A. (24) 

In this case it is well known that the usual Darboux scheme works and that the Fj are 
actually A-independent functions Fj{\,x) = fj{x). More precisely, as explained above, we 
must choose particular eigenfunctions '■Pjivj) of the Schrodinger operator d1 + Wj{x) (with 
eigenvalues Vj) such that the fj{x) are expressed as 

fj{x) := (log(^j)^, (25) 

hence these latter functions have to satisfy the following specialization of relation 

fi{x)' + f^{xf + w,{x)-u, = Q. (26) 

Simultaneously we have to satisfy equation (jH)) and therefore we choose without loss of 
generality, 

/i,(A) = (z/,-A), A,{x) = l. (27) 

This then yields the standard Darboux transformation for the Schrodinger operator with 
fj = fj (from ©), Gj{x) = [dx — fj{x)) (from (j2I)), and with a change in the potentials 
dZj) given by: 

Wj+i = Wj + 2fj = Wj + 2 (log (fj)^^ . (28) 

The resulting chain equation ()10|) for such Darboux transformations is the well known 
dressing chain [T] 

fj + f'j+i = ff - f]+i + aj, aj = T^j+i - yj- (29) 

The linear problem for the dressing chain [2Sj follows from the general expressions 
(fT2lT3|l and (da in, subject to (jSIj): 

[dx-S-fj]^, = (30) 

[S^ + f,S + Sfj + v,-\\ij, = 0. (31) 

From equation (fTBj) it can be seen that the compatibility condition of this linear system 
takes the form of a straightforward operator commutation relation [M^, L^] = 0. 
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3.2 Generic bilinearization 



Note that the initial potential, say wq, in the sequence never appears in the dressing 
chains. In fact, it only becomes relevant when actually (re-) interpreting the solutions of 
these chains - through - as solutions of Schrodinger equations such as (pUj) . which 
explicitly depend on specific potentials. One instance where the underlying Schrodinger 
equations are of importance is when one wishes to obtain a Hirota bilinear form for the 
dressing chain. For, if in the sequence of coupled Schrodinger equations 

(dl + Wj)ipj = Ujcpj, {] = 1, 2, 3), (32) 

with potentials transforming as in ()28j) . we parameterize each potential Wj as 

Wj = 2(logu;j_i)2x, 

then we obtain a multiplicative transformation rule for the new functions Uj 

This suggests a parameterization of the eigenfunctions ipj (j32|) as the ratio of uj and ooj-i. 
As is well known, the standard bilinearization of such Schrodinger equations is through a 
ratio of r-functions. Here, with the benefit of (considerable) hindsight we set 

and then (jH^ transforms into 

{Dl - exD^ - kj^i)t, ■ Tj^i = 0, (34) 

for some constants e and kj. (For a definition of the Hirota /^-operators and an introduction 
to their importance in the context of the Painleve equations we refer to the review paper 

It can be shown that this chain of bilinear equations is nothing but a similarity reduction 
of the (2 + l)-dimensional dressing chain associated to the modified KP hierarchy. This 
entitles us to refer to these Tj as genuine tau-functions in the sense of Sato theory [221123121] ; 
the specific description of this reduction will be addressed in a separate publication pUj . 
Note that, as pointed out above, the initial tau-function tq (and hence the potential Wq) 
appears explicitly in the bilinear chain (j34|) . 

Below we will see that the system (jH^ provides a generic bilinearization for all Painleve 
equations (and higher order variants) contained in the dressing chain. There do exist other 
possibilities when it comes to bilinearizing the Painleve equations, see [23 I2S1 1101; but 
we believe the present generic approach to be new. (For a general survey of the Painleve 
equations as similarity or symmetry reductions of integrable nonlinear partial differential 
equations, see, e.g., PT|.) 

3.3 Periodic closing of the dressing chain 

As was mentioned in the introduction we are interested in dressing chains of finite length. 
Specifically, we will impose periodic closing conditions on the sequences of Darboux trans- 
formations used in the construction of the chains, i.e., following Adler jB], we require that 
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- up to a shift in the eigenvalues Vj - the eigenf unctions <^j{i^j) which define the Darboux 
transformations (as in ()25|) ) become periodic with some period > 1: 

iPj+N{i^j+N) = with z/j+Af = z/j - e, (35) 

for a shift e ^ 0. 

As these functions satisfy {d^ + Wj)ipj = I'j^j, we immediately find that the periodic 
closing of the sequence of functions ipj also implies that the sequence of potentials Wj 
(generated through these very Darboux transformations) has to close as 

Wjj^N^x) = Wj{x) — e, (36) 

with the obvious implication that the sequence of generalized potentials Uj{X, x) closes 
as 

Mj+Ar(A, x) = Uj{X + e,x), (37) 

for generic eigenvalues A. 

The closing conditions (j35|) take on a particularly simple form when expressed on the 
bilinear chain (jHH): 

Tj+N = Tj, Kj+N = i^j, and e = (38) 
It should be remarked that as a consequence one also obtains the constraint 

N 

\{^,=e-^^''\ (39) 
i=i 

from (jHHj) . Finally, for the variables (and parameters) which appear in the dressing chains 
()29|) . the periodicity conditions take the form: 

N 

fj+N = fj, aj+N = aj, ^aj = -e, (40) 
accompanied by the first integral (which holds for all > 0) 

N 

If A^ is even there is another integral, obtained by summing every other term of ()29p : 

N N 

=-^E(-i)'«^- (42) 

Note that the cases A^ = 1, 2 are therefore solvable by quadratures. 



4 Piv 

As was discovered almost a decade ago j21 E] , the 3-periodic dressing chain (pUj) is nothing 
but the Pjv equation. We shall now go on to show this explicitly on the dressing chain, 
after which we shall derive a bilinear representation and a Lax pair for the P/y-equation. 
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4.1 The symmetric form of Pjy 

Using condition (j4(J|) at = 3 we obtain from (j29j) the chain of equations 

/{ + /2 = /f-/l + «l, 

/2 + /3 = /|-/| + «2, (43) 
+ =/|-/f + «3, 



with 

Expressed in terms of 
we get 



«! + ^2 + tts = (44) 

91 = fl + f2 = (logV?lV22)x, 

92 = f2 + h = (logv?2V23)x, (45) 

93 = f3 + fl = (logV33V3l)x, 

9'i = 9ii93 - 92} + tti, 

^2 = ^2(^1 -^3) + "2, (46) 
^3 = ^3(^2 - ^1) + as- 
Following Noumi and Yamada fSJ QHl we shall refer to this system as the symmetric form 
of the Pjv equation. This particular form of the period 3 dressing chain appears already in 
j3] . However, it is worth pointing out that (pUj) is already presented in where it is used 
to integrate a higher order non linear differential equation in terms of the Piv equation. 
Because of the constraint - or alternatively, as a consequence of pij) - the system 
can be integrated once: 

91+92 + 93 = -ex (47) 

(where a possible integration constant only amounts to a translation in x and can therefore 
be omitted). Eliminating we get 

g[ =ai- 2gig2 - exgi - gf, 
g'^ = a2 + 2gig2 + exg2 + gl 

These equations are in Hamiltonian form: 

/ , OH 

9i = 7^, 92 = -j^, 49 
dg2 dgi 

where 

H = -gigl - g2gl - exgig2 - a2gi + aig2. (50) 

Note that this Hamiltonian differs slightly (by a canonical transformation) from the one 
discussed by Okamoto [3 HE] ; it does, however, appear in this form in the classic work 
inS], and in dS]. 

The Piv equation is now obtained if we eliminate g2 and denote y{z) = K,gi{x), where 
2/e, this yields the standard form 

(fy I fdyY 3 . 2 , Of -2 ^ , ^ ^i^i^ 
+ -y +4:zy +2{z - a)y + - (51) 



dz^ 2y \dz J 2 y 
where the parameters a and b are given by 

a = {a2 — a3,)/e, h=—2{ai/eY. (52) 
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4.2 Bilinear form of P/y 



Imposing conditions (j38j) on the bilinear form (j34j) at = 3 we immediately obtain a 
bilinearization of the P/y equation: 

(Dl - fD^ - Ko)ri ■ To = 0, 

{Dl - fD, - Ki)r2 ■ n = 0, (53) 
{Dl - fD, - K2)to ■t, = 0. 

The transformation for the Qi is: 

Qi = log ^ - ai = K,i- Ki_i - J, (54) 

Ti^i 3 3 

for i = 1,2, 3, and with periodicity (fH^jl (see also 15, 16j). Recall that was obtained by 
making the substitution (jH^ directly into the (now periodic) chain of Schrodinger equations 
()32j). If, on the other hand, we make the substitutions ()54|) into the symmetric form ()46p 
we will only obtain two equations for the three tau-functions. This is because the equations 
do not contain any information on the potentials w. However, if we explicitly include 
the assumption Wj = 2(logrj_i)2x, e.g., in the form of the first equation of (jSSl), then we 
do obtain a well-determined system of bilinear equations. Note also that 

y{z) = d Jog (^^e-'r^ . (55) 



4.3 Lax pair for Pjy 

Let us see what the consequences of the periodic closing (j35|) are on the linear problem 
()30I31|) . According to the generalized potentials Uj{\,x) are periodic for generic A as 
well (up to a shift in A) and thus we may require similar periodicity for some eigenfunctions 
of the linear problem (^, namely those generated by successive Darboux transformations: 

tpj+3{X,x) = i!j{X + e,x). (56) 

The linear system ()30|31|) therefore reduces to the following set of "difference equations" 
in the spectral parameter A: 

/ipi\ fi>2 + fiipi\ f i'3 + 9ii'2 + M\ 

5. ^2 = U3 + /2^2 , 92^3 + ^2^2 ] = X { ^2 ] , (57) 

VW VA' + ZsW \r2+93ri + 'y3^3j vW 

where the symbol tpj denotes eigenfunctions at shifted values of the spectral parameter: 
iljj{X,x) := i/jj{X + e,x). 

In order to recover a more customary linear problem for the Pjv equation, let us first 
rewrite the equations (fHTjl as: 

d^'^{X) = Bi^{X) + B2^{X + e), (58) 
(A-AI)^(A)+^2^(A + e) =0, (59) 
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for ^(A) = {ipi{X, x) , ■iIj2{^, x) , ips^X, x)y and with matrices 81,82, Ai and A2 



ffi 1 0\ /O 0\ 1\ fO 0\ 

81= I f2 1 , i32 = , ^1 = U2 92] , A2 = { 1 0. (60) 
\0 /g/ Vl 0/ yO u,J \g, 1 0/ 

If we now define the (formal) Fourier transform of ^I/(A) as 

^(k) = [ d\ e'^^ ^(A), (61) 



we immediately see that $(A;) satisfies local relations in k, for example: 

a^l'(fc) = 8i^{k) + i32e-^='=l'(A;) (62) 



instead of the original difference relation ()58|) in A. Consequently, introducing a new 
parameter 1^ and a new function $(0 in terms of $(/c) by 

^=e-^ $(0:=<^'(^) . , (63) 

we find that $(0 satisfies the following Lax pair for the Pjv equation, 

d^^ = M<^, M = 81+^82, , 
-eid^^ = £$, C = Ai + iA2, ^ ' 

and the compatibility condition dxC + eE,d^M. = [Ai,C]- is nothing but the symmetric 
form PU)) for Piy- This Lax pair already appears in jTH], in connection with reductions of 
the Drinfeld-Sokolov hierarchy. Also, the transformations ()6H1 and ()6H|1 combined, amount 
to the Mellin transformation that is used (to the same effect) in j21j in order to connect 
Lax pairs arising in the context of dressing chains to those that appear in the work by 
Noumi et al. 



5 Pv 

Similarly to the above, Py can be obtained by closing the dressing chain ()29p at = 4 
13 El • Just as for the Pjy equation this insight yields immediate access to a Lax pair and 
bilinear formulation for the Py equation. 

5.1 Periodic closing and Py 

Closing the Darboux chain (f^^ at period = 4 (cf. condition ()40|) ) yields the following 
system of differential equations: 

/{ + /2 = /f-/l + «l, 

/^ + /H/|-/| + «2, .... 
/^ + /l = /|-/l + «3, ^^^^ 
+ = /|-/f + «4, 
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with constants aj {j = 1,2,3,4) subject to the constraint 

«! + ^2 + «3 + ct4 = (66) 

Now we have two first integrals 

/1 + /2 + /3 + /4 = -lex, (67) 
-/1+/2-/3+/4 = i(«i-«2 + «3-«4)=:c^, (68) 

(recall ()4ip42j) ). Hence, the system (jHSI) can be reduced to a 2'^'^ order differential equation, 
which turns out to be nothing but the Painleve V equation. We shall now proceed to show 
this. 

It is convenient to rewrite the symmetric form ()65|) in terms of new variables gi,g2 
obtained by resolving (|U7j) and (|UHj) as 

/s = -^9i92- ligi- g2) - f + ^, 
fA = i^9i92-\{9i + g2)-f-^, 
the inverse relation of which is 

gi = /l + /2 - /3 - /4, 

g2 = -/i + /2 + /3 - h, 

together with the constraints (jHTjl and (jHH|l . In terms of these new free functions gi and g2 
we obtain the symmetric form 

g[ = ~9l92 + %9i + "-f92 + a^-a^, (69) 

92 = i:gigl-j^g2-fgi + a2-a^, (70) 

which is Hamiltonian with 

H ■= -^9l9l + f{9l + 9l) + 7^9192 - («2 - ai)gi + ("i - a^)g2, (71) 

Note that this Hamiltonian differs from the form used in [23 13 IQ but they are 
connected by a canonical transformation and a change of the independent variable. 

The Py equation is now obtained if we first solve g2 from (jU^ . substitute it into (fTOI) 
and express gi in terms of y defined by 

y = \~ 7^91^ where y = y{z), z = ex^/A, 

this yields the Py equation in the form (used, e.g., in [221122], up to an extra transformation 
of the independent variable) 

(fy f I ^ 1 ^ f dy\^ Idy 



where 



dz'^ \2y 2{y — 1) J \dz J z dz 

_ ay _ l3{y - 1) _ -/yjy - 1) _ 
z'^{y — 1) z^y z 

(ai — a-i — uj + ie)^ (ai — + u — ^e)'^ 



Syiy-l){2y-l), (72) 



otj n [Oi - 03 + uj - ^8)- ai-a2 . i 

The usual form of Py is obtained for w = y/{y — 1), which permutes the poles at 1 and 00. 
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5.2 Bilinear form of Py 



Imposing conditions (j38p on the bilinear form (j34j) of the dressing chain, this time at = 4, 
we obtain a bihnearization of the Py equation in terms of 4 tau-functions: 

{Dl - f - Ko)ri ■ ro = 
{Dl - fD, - k,)t2 ■ n = 

{Dl - fD, - K3)to ■t, = 

Subject to the same remarks as in the case of the Pjv equation, this system of Hirota 
equations gives a bihnear representation for the symmetric form using substitution 

fi = d, (^og - y , ai = Ki- Ki_i - |, (74) 

for i = 1,2,3,4, and with periodicity Note also that (cf. (fS3j)) 

y{z) = d Jog (^e^^ . (75) 

It was already mentioned that in ^3] Noumi and Yamada presented a class of dynamical 
systems, each member of which possesses a particular symmetry group of (affine) Weyl- 
type W{A^^}_^). Pjv and Py are contained in this class at levels n = 3 and 4 respectively. 
In can be shown [331130]! not only that the bilinear forms of the P/y and Py equations f (j33jl 
and ()73|) ) possess similar symmetries, but generally that any bilinear system contained in 
(|3^ possesses a symmetry group of type W{A^j^}_-^). Hence these equations are nothing but 
the bilinear formulations of the Noumi- Yamada systems and consequently, we have shown 
that each member in that class corresponds to a periodic dressing chain (in the sense of 
Sects. 3.2 and 3.3). 



5.3 Lax pair for Py 

Exactly as in the P/y case one can obtain a Lax pair for Py by imposing periodicity on 



the generic Lax pair 



for the dressing chain: 



Following the notation of Sec. 14.31 we can write: 



fi^3 + (/l + f2)^2 + ^^l^A 
^4 + if 2 + fzj^l. + V2lp2 

V'l + (/a + /4)^4 + v^ipz 

\r2 + (/4 + hWl + ^^4^4/ 



^2 
V^4/ 



After a (formal) Fourier transformation similar to the one in the P/y case 
obtain the following Lax formulation for the Py equation |ltij : 



(76) 



(77) 



, we 



(75 



13 







ffi 1 








/o 


o\ 










M = 


/2 1 























/s 


1 


+ 4 


















^0 


fj 








0/ 










/l + /2 1 





\ 






/ 








o\ 







/2 + /3 


1 




+ e 






















/3 + /4 


1 



















Z/4 


J 
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(79) 



^0) 



The compatibility condition d^C + e^d^M. = [A4,C]^ is easily seen to correspond to the 
symmetric form of the Py equation. 



6 The Pjii branch in general 

As was noted in the introduction, the Pm branch is quite different from that for the 
Pjv — Py equations. In this section we will first study this branch in general and in the 
next section we shall close the chain in order to obtain Pni- 

6.1 The generic chain equations 

As was explained in Section IT^ a different branch of the chain equation ()10|) for potentials 
of type (fT7|) is obtained if we choose 



(81) 
(82) 



F,{\x) := \ + fj{x), 
/i,(A) := 2(A-z/,), 

A,{x) := [f^{x) + \v^{x)]-\ (83) 
(/j + itij ^ 0); the only remaining part of condition (jH)) is the following constraint on / 



(84) 



For the present transformation to be a genuine Darboux transformation, we must also 
require that the intertwiner Gj in Q annihilates some particular eigenfunction ^pj : = 
ipiyj^x) of the original Schrodinger equation, i.e., Fj{h'j,x) = (\ogipj)x or in terms of /: 

fj ■= -^j (85) 
Condition ()84|) then turns into the Schrodinger equation for the eigenfunction (pj: 

V"i + (-'^i + + = 0- (86) 

The transformation induced for the potentials is obtained from (|7j): 



Vj+i = Vj + 2{\og Aj)' , 

wj+i = w, + [2fjA,-A'.]'/Aj. 



17) 
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From the explicit form (j83|) of A it is clear that in the present case the chain equations 
generated by ()10|) - subject to condition (jH^ - will take the form of a dynamical system in 
the variables / and v. However, these variables are strongly coupled and from (j84l87l88|l 
we get (after eliminating w) 

f'j + ¥3 + + - ^i) = 0, (89) 

/i - H - (A-i + '^.-i)' + (¥-1 - + (/, + ^jf - {\v, - = 0. (90) 

It is convenient to change to new dependent variables dj, rj defined by 

f3 = ¥j + ¥j-^J' "^j = dj-rj+2uj, = z/j - z/j+i, (91) 

(note that dj = fj + ^Vj = Aj"^) because then the chain equations of the Pm branch take 
the following simple form (which already appears in ^) 

dj = dj{dj - rj - dj+i + rj+i + 2(3j), (92) 
r/ = dj.irj^i- djVj. (93) 

The linear problem for this integrable system is of course obtained from p2pi3|l and 
()15I14|) by appropriate identifications of the variables appearing in those equations, but we 
shall not discuss its general form here. 

6.2 Generic bilinearization 

It is of interest to bilinearize the chain equations ()92f93|l before closing the chain. We start 
by representing A in terms of a new function (fj and the eigenfunction (pj-. 

Af = ^, (94) 
^j 

which due to (|HHjl amounts to 

= (pj + {vj - \vj)ipj. (95) 
Since ipj is a solution of (jHBj) we also obtain 

Taken together, equations ()95p9(ij) yield a first order representation of the Schrodinger 
equation (jH!)|l and as such provide a suitable starting point for the bilinearization of the 
Schrodinger equations associated with the Pm branch. 

Under the above ansatz, the transformation formula ()87p suggests a logarithmic pa- 
rameterization for the potential v 

vr.= {\ogq,)\ (97) 
because then the transformation for q takes the simple multiplicative form 

~ (98) 
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where ~ denotes equivalence up to a constant multiple. The transformation formula for w 
is best expressed in terms of a new variable W 

= l^j - (99) 

because then we get 

= w, + (log - (log q,y + (log (p,y\ (loo) 

(from (jHHj) using also ()95I98|) ) which suggests a parameterization 

W, ■= (logo.,)". (101) 
The transformation for W then takes the simple multiplicative form 

u,^, ^ ^ 5^ (102) 

(where ~ denotes equivalence up to a gauge factor e"^"'"''). 

If, as before, we think of the functions ipj, (fj, qj and Uj as ratios of "tau-functions" , it 
becomes clear that 

= e^^^, = T,T^, (103) 



'j'j ' 




cpj = c,^J-^ e"^", ipj = g . ^ W ^ e"^", (104) 

V V 

is a good representation, since it automatically satisfies the multiplicative transformation 
rules for v and equations and ()102|) . if the introduced constants satisfy 

7i+i = 7i + ai -ttj- (105) 
For / and v we get from (|H^ and (jHTj) 

/, = log (^-^^ + - z/,, = log (i) + 7i. (106) 

and then 

dj = log (^^) + + l-fj - Uj, Tj = log + aj - + z/^-. (107) 

As mentioned before, the assignments ()103I104|) take care of equations ()87I88|) . but we 
still have to study the Schrodinger equation (IHU]). which we take in its first order form 
dnSEEl). This yields 

(/},. + a, - V, + f ) r/^i ■ r, = I e^^^-^^^^ r/r.+i, 

D.(D. + a, - V, + ^) r/^i • r, = I e(^^-<^^)- (D. + 7, - - - 2z/,)r; ■ r,+i. 

In p03ll04|) we have some additional freedom in the parameters 7, a, a (by scaling r^^ ^ 
^+gSjx^ Tj — > Tje"**^"^'), which can be used to assign a fixed value to the parameter 7^ 

7, = 2(z.,-a,). (108) 
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In order to satisfy p05|) we must then also require that 

(kj = 2aj+i — aj — 2(z/j+i — Uj). (109) 
After this redefinition of parameters, the equations simphfy considerably and we get 



where 



Pj = "i + 3«i = 2(aj+i + aj - Uj+i + z/^), 



K 



Cj /Cj . 



Equations piO|) form the Hirota bilinear representation of the chain equations ()92|93|) . For 
later use note that due to the first equation of piO|) we can also write 

d, = e'^^-i^. (Ill) 
6.3 Periodic closing in general 

As before, we are mainly interested in a finite chain of equations obtained after requiring 
the potential Uj{\) to be periodic, up to a shift in the (generic) eigenvalue, as stated in 
(|?fj). and this implies the existence of "periodic" eigenfunctions For the potentials 

V, w this now implies 

Vj+N = Vj — 2e, Wj+N = + evj — e'^, (112) 

whereas the auxiliary potential W of is strictly periodic: Wj+n = Wj. 

On the level of the r-functions r, r"*" these closing conditions translate into strict peri- 
odicity 

r^+N = T-/, rj+N = Tj. (113) 

Condition p()5|l was replaced by constraints p()8|l and p()9|l . which are compatible with 
the requirement that the parameters aj, aj, cj, Cj, aj, pj, kj are strictly periodic and 

jj+N = lo - 2e. (114) 

We also have 

N N 

^/5,=£, 5^ a, = 25. (115) 

Furthermore, from ()107ll08p it is obvious that dj and are (strictly) periodic as well and 
from ()111|) and ()107p we also find the following two (generic) conservation laws for the 
periodic chain equations jl]: 



N N N N 

n d, = e^- n E = 2 E (116) 

j=l j=l j=l j=l 
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7 Pn, 



7.1 Closing the chain at = 2 

Due to the existence of two conservation laws the system obtained from (|92l93p at = 1 is 
of course trivial. However, at = 2 we find that di,d2,ri, r2 satisfy the 4 basic equations 

di' = di{di-d2 + r2-n + 2pi), (117) 

d2 = d2{d2-di+n-r2 + 2(32), (118) 

r/ = (^2^2 - diTi, (119) 

r2' = diTi - d2r2, (120) 

with Pi + f32 = £■ We shall refer to this system as the "symmetric form of Pm" . We now 
show how Pjjj is obtained from it. 

Two first integrals were given in pi6|) 

did2 = KiK2e^^'', n + r2 = |(pi + P2) -e (121) 

We resolve these by introducing new functions Qi defined as follows: 

di = Ki gi e^"", d2 = 1^2 Qi^ e^""', n = -gig2 + \{pi + P2) ^"2 = 9192- (122) 

Then the equations ()117lll20|l take the form 

g[ = 2glg2 + e'^mgl-{p2-e)gi-e'^K2, (123) 
g'^ = -2gigl-2e''-'Kigig2 + {p2-e)g2 + e'''Ki{l{pi+P2)-e), (124) 

which is Hamiltonian with 

H = gfgl + e'^'Kiglgi - (p2 - e)^i^2 - e^^[/ti(i(pi + P2) - e)gi + ^2^2] (125) 

(which is the same as the one given in [HEIIEI; up to a simple transformation). 
If we solve g2 from ()123p . substitute it into ()124p . and use the new variables 



y{z) = Jf^gi{x), where z = V^e^, (126) 



we get Piii in the canonical form: 



dz^ y^dz^ zdz^^^z^^^ y ^ ^ 

7.2 BiUnear form of Pm 

The bilinear form of Pm is basically nothing else but ()110j) with periodic r's. Since the 
r's appear in pairs (ri, and r2, t^) and since bilinear equations are gauge invariant, we 
propose the following gauge transformation and scaling 

ri = v^rie''^ rt = fte'\ r+ = f+, t<, = ^^f^, ^ = i(^^ _ a^), (128) 
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which yield, after a subsequent x-translation x — > x — ^ log(fi;i/t2), the following novel 
bilinear form of Pi if. 



D,f+-n = e^^f+f2, (129) 

D.,ft-f2 = e'^f+n, (130) 

Dlf^-n = e^^(D,-pi)f+-f2, (131) 

Dlft-f2 = e'"" (D^ - p2)f+ ■ n. (132) 



From this Pm can also be derived directly: Comparing ()122I111I126|) (after the aforemen- 
tioned gauge transformation) one can solve for f2 in terms of y and the other f's. Then 
from (jl29|) one solves for d^fi and from (jl3(jp for d^f^- Then (jl31|) x^— (jl32j) is nothing 
but the Pi 1 1 equation ()127|) . In these new variables y and z we also have (cf. ()55I75|) ) 

yiz)=dAog^. (133) 

n 

It should be noted that bilinear forms are sensitive to changes in the independent 
variable. For example, if we express the bilinear system obtained after the transformation 
()128|) (but without the x-translation), in terms of the z-coordinate (now regarding the tau- 
functions f2 ,Ti,fi ,T2 as functions of z) we get a system of equations which is no longer 
expressible in Hirota D-operators only as it also involves ordinary z-derivatives. However, 
the system so obtained can be shown to be a reduction of an integrable system contained in 
the so called "modified 2-component KP hierarchy" (at least in the restricted case where 
P2 = ±Pi) and hence the functions f2 ,Ti,ft ,T2 introduced here can be thought of as 
genuine tau-functions in the sense of Sato theory. 

7.3 Lax pair for Pni 

Just as for Piv, a Lax pair for the Painleve III equation can be derived from the linear 
formulation p3ll4|) introduced earlier for the general chain equations (fTUI) . 

As before, the standard periodic closing (j35|) involves a shift in the eigenvalues Uj 
associated with the eigenf unctions ipj (for the Pm case one has = ui — e). Hence, from 
the condition ip3{i'3) = and expression (jH3j) it follows that /s = fi+e. Furthermore, 

due to the relation u^^X) = ui{\ + e) there will exist eigenfunctions of the Schrodinger 
equation for which ■?/'3(A,x) = ipii^ + ^j^) holds. This then, bearing in mind the 
relation (|HT|) and the fact that = Ai (actually A^ = Af ; we choose the sign A^ = Ai), 
yields the following periodic closing of the linear equations ffT^ 

9,^1 = ArV2 + (A + /i)^i (134) 

^,^P2 = A-Vi + (A + /2)V'2, (135) 

and (HH) 

(AiA2)-Vi + ^r'(2A + /i + /2-(logAi)')V'2 + 2(A-z/i)A^2^i = 0, (136) 

(AiA2)-V2 + ^2'(2A + £ + /i + /2-(logA2)')V^i + 2(A-z/2)A~'^2 = 0, (137) 

where ip^ stands foYipi{X + e,x). 
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It is now advantageous to change to scaled eigenfunctions 

(recall that ri + r2 is a constant). After also changing to the new independent variable 
^ = ez = y/KiK^e^^ and using the previously obtained formulae we can write (jl34H137j) as 

e^a^^i = i(4A + rfi + rf2 + e-4z/i+4z/2 + 2ri-2r2)^i + V3^2, (138) 
£3^3^2 = i(4A + c/i + c?2-£ + 4z/i-4z/2-2ri + 2r2)z^2 + v^V'i, (139) 



iijl + y^{2\ + di + r2 - 2u2)i^2 + 2di(X - = 0, 

3V'2 + V3(2A + c/2 + ri -2z/i + 2£)^^ + 2d2(A-z/2)V^2 = 0. 



(140) 
(141) 



In exactly the same way as for the Pjv equation, this system of difference equations (in 
the spectral parameter) can be cast into a more standard form by introducing the formal 
Fourier transform of the eigenfunctions il!j{X,x): 



dX e'^^ 



ijj{x,x) 



(142) 



In terms of the new variable ^ = exp{—iek) and the new dependent variables $^(^,3) : = 
^j{k,x), we then obtain the linear systems 



$1 
$2 



^3^3 ( 



Di + A + 
D2 + A- 



1 



1 



$1 

$1 
$2 



where 



A 



-idiTi/^ ir2 ) 



(143) 
(144) 

(145) 



D, = -(z/i + z/2)( J ? ), D2 = [lid^ + d2) + u, + U2]{h^,) + [^s + n-r2]a -°i )• (146) 
The compatibility of these two matrix equations yields the equations pi7lll2()(l . 



8 Conclusions 

In this paper we have constructed Darboux chains from (scalar) Schrodinger equations for 
generic second order energy-dependent potentials and then proceeded to discuss periodic 
reductions of such chains. The Darboux chains were classified into three types, the so called 
Piv-v, Pi 1 1 and Pyi branches, which derive their names from the Painleve equations that 
arise as the lowest period (non-trivial) reductions contained in each of them. (A detailed 
discussion of the Pyi branch will be given in a sequel to the present paper.) 

We described in detail the construction of the generic chain equations (j29|) and ()92p93|) 
for the Piv-v and Pm branches, and their bilinearization ()34|11()|1 . For the periodically 
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closed chains we described the reduction to the corresponding Painleve equations, with an 
intermediate equation in Hamiltonian form (see also |3]). 

The bilinearization of the generic Darboux chains in the Piv-v and Pm branches 
led to a representation of the eigenfunctions in terms of tau-functions (j33ll(J4j) . which in 
reduction gave rise to bilinear representations and tau- function formulae for the Pm-v 
equations ()53I73I129II132|) and their solutions ()55I75I133|) . From the bilinear form of the 
Piv-v branch it is clear that the equations described by it are identical to the so called 
An ''-type dynamical systems introduced by Noumi and Yamada [14. The exact nature 
of the tau-functions associated to the periodic reductions of this branch in the context of 
Sato theory will be discussed in a forthcoming paper (see [M] for a discussion of the Pjv 
case). The precise link between the tau-functions that appear in the periodic reductions 
of the Pjii branch and Sato theory is currently being investigated. 

The construction of the Darboux chains presented in this paper also allowed us to 
systematically derive Lax representations for the chain equations and their reductions, 
resulting in a novel Lax description of the Pni equation ()143II146|) . The Lax pairs obtained 
for the Piv-v branch were again related to those obtained by Noumi and Yamada for the 
An^^-type dynamical systems mentioned before. 

The interpretation of the Hamiltonians associated with the periodic reductions of the 
Darboux chains, in terms of the tau-functions that describe their solutions, as well as a 
detailed investigation of the higher order members - i.e., those corresponding to chains 
with higher periods - in the different branches of chain equations, are topics that will be 
addressed in subsequent papers. 
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